Theory of proximity effect in superconductor/ferromagnet heterostructures 



A. Bagrets and C. Lacroix 
Laboratoire Louis Neel, CNRS, B.P. 166, 38042 Grenoble cedex 9, France 

A. Vedyayev 

CEA/ Grenoble, Departement de Recherche Fondamentale sur la Matiere Condensee, SP2M/NM, 38054 Grenoble, France and 
Department of Physics, M. V. Lomonosov Moscow State University, 119899 Moscow, Russia 

(Dated: February 1, 2008) 

We present a microscopic theory of proximity effect in the ferromagnet /superconductor- 
/ferromagnet (F/S/F) nanostructures where S is .s-wave low-T c superconductor and F's are layers 
of 3d transition ferromagnetic metal. Our approach is based on the direct analytical solution of 
Gor'kov equations for the normal and anomalous Green's functions together with a self-consistent 
evaluation of the superconducting order parameter. We take into account the elastic spin-conserving 
scattering of the electrons assuming s-wave scattering in the S layer and s-d scattering in the F lay- 
ers. In accordance with the previous quasiclassical theories, we found that due to exchange field 
in the ferromagnet the anomalous Green's function F(z) exhibits the damping oscillations in the 
F-layer as a function of distance z from the S/F interface. In the given model a half of period of 
oscillations is determined by the length ^ n — tvvf/sc*, where vf is the Fermi velocity and e cx is the 
exchange field, while damping is governed by the length lo — (l/l] + with 2-j and l± being 

spin-dependent mean free paths in the ferromagnet. The superconducting transition temperature 
T c (c(f) of the F/S/F trilayer shows the damping oscillations as a function of the F-layer thickness 
dF with period £f = ir/y?me ax , where m is the effective electron mass. The oscillations of T c (c(f) 
are a consequence of the oscillatory behavior of the superconducting order parameter at the S/F 
interface vs thickness dF, that in turn is caused by the oscillations of F(z) in the F-region. We show 
that strong spin-conserving scattering either in the superconductor or in the ferromagnet signifi- 
cantly suppresses these oscillations. The calculated T c (<1f) dependences are compared with existing 
experimental data for Fe/Nb/Fe trilayers and Nb/Co multilayers. 

PACS numbers: 74.80.Dm, 74.50.+r, 74.62.Bf, 74.76.-w 



I. INTRODUCTION 

The artificially fabricated layered nanostuctures with 
alternating superconducting (S) and ferromagnetic (F) 
layers provide a possibility to study the physical phe- 
nomena arising due to proximity of two materials (S and 
F) with two antagonistic types of long range ordering. 
One of such interesting effects is the existence of so- 
called 7r-phase superconducting state in which the order 
parameter in adjacent S-layers has opposite sings. The 
7r-junctions were originally predicted to be possible due 
to spin-flip processes in magnetic layered structures con- 
taining paramagnetic impurities in the barrier between 
S layers^. Later on, Buzdin et alm£ and Radovic et alX 
showed that, due to oscillatory behavior of the Cooper 
pair wave function in the ferromagnet, 7r-coupling can be 
realized also for S/F multilayers. The 7r-coupling leads 
to a nonmonotonic oscillatory dependence of the super- 
conducting transition temperature T c as a function of 
ferromagnetic layer thickness dp 2,3,4 . The effect occurs 
because of periodically switching of the ground state be- 
tween 0- and 7r-phase, so that the system chooses the 
state with higher transition temperature T c . 

These theoretical predictions stimulated a considerable 
interest to proximity effect in S/F structures also from 
experimental point of view. First the oscillatory behav- 
ior of T c (cIf) was observed by Wong et alA in V/Fe mul- 
tilayers and later on these results were well explained 



by theoretical calculations of Radovic et alA However, 
in subsequent experiments with V/Fe multilayers 6 , the 
oscillatory T c (dp) dependence was not observed. The 
following experimentsi^iiSiiiiiSiH revealed the different 
and even controversial behavior of T c {oIf) for different 
structures. The nonmonotonic oscillation- like behavior 
of T c (cIf) was reported by Jiang et al£ for Nb/Gd mul- 
tilayers and Nb/Gd/Nb trilayers, by Miihge et al& for 
Fe/Nb/Fe trilayers and recently by Obi et al£ for Nb/Co 
and V/Co multilayers. However, negative results were 
reported for Nb/Gd/Nb trilayers by Strunk et al*&, for 
V/Vi-zFej; multilayers by Aarts et al. 11 , for Fe/Nb bi- 
layers by Miihge et alw£ and Nb/Fe multilayers by Ver- 
banck et alw&. For interpretation of experimental results, 
along with mechanisms of 7r-coupling and suppression of 
T c due to strong exchange field in the ferromagnet, an- 
other mechanisms were suggested such as complex be- 
havior of the "magnetically dead" interfacial S/F layer 
(see details in Ref. @), the effects of a finite interface 
transparency^, and spin-orbit scattering 14 . 

The original theory of proximity effect proposed by 
Buzdin et alm£ is based on the quasiclassical Usadel 
equations^ applied for S/F structures. In this case the 
Usadel equations must be supplemented by boundary 
conditions for the quasiclassical Green's functions at the 
S/F interface. This essential point was recently discussed 
in Ref. ^| On the other hand, the boundary condi- 
tions for microscopic Green's functions can be written 
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obviously for ideal S/F interfaces if one uses Gor'kov 
equations 1 ^. These equations, however, are more com- 
plex to resolve than the quasiclassical ones. In the given 
paper, we present a theoretical investigation of T c (d F ) 
behavior for F/S/F trilayer structures based on Gor'kov 
equations. We consider that F layers are 3c? transition 
metals and assume that the main mechanism of spin- 
conserving electron scattering in F layers is s-d scatter- 
ing, while S layer is s-wave superconductor with s-s scat- 
tering. We find the characteristic lengths determining the 
periods of oscillations and damping of critical tempera- 
ture T c and Cooper pair wave function, and show that 
in the given model these lengths differ from length scales 
predicted by quasiclassical theorie a 3 ' 4 ^' 1 ^ . We show that 
strong spin-conserving scattering either in the supercon- 
ductor or in the ferromagnet significantly suppresses the 
oscillations of T c . We compare our results with exist- 
ing data on T c (d F ) for Fe/Nb/Fe trilayers& and V/Co 
multilayers^, where F's are 3d ferromagnets, and find rea- 
sonable agreement with theory and experiment. 



II. GOR'KOV EQUATIONS AND GREEN'S 
FUNCTIONS 



We consider a trilayer structure F1/S/F3 where S is 
low-T c superconductor and F's are 3c?- metal ferromag- 
netic layers. The thicknesses ds and d F of the S- and F- 
layers are supposed to be much smaller than the in-plane 
dimension of the structure, so that the system can be con- 
sidered as homogeneous in the xy-plane (parallel to the 
interfaces). We denote the axis perpendicular to the xy- 
plane as z-axis. Let z = ±a be the positions of the outer 
boundaries of the F-layers and z = ±d be the positions of 
S/F interfaces, then c?s = 2c? and d F = a — d. We adopt 
that S is a simple s-wave superconductor with s-s mecha- 
nism of electron scattering. According to Ref . I22I for su- 
perconducting Nb which is usually used in preparing the 
S/F heterostructures, s-wave scattering is indeed prevail- 
ing. Concerning the ferromagnetic layers, we adopt the 
simplified model 23 considering that two types of electrons 
form the total band structure of 3c? transition metals: al- 
most free-like spin-up and spin-down electrons from sp 
bands (these electrons are referred as s electrons) and 
localized d electrons from narrow strongly exchange split 
bands. The main mechanism of spin-conserving electron 
scattering in 3g? ferromagnetic metals is s-d scattering 24 
because of a dominant contribution of d density of states 
(DOS) to the total DOS at the Fermi energy e F . The 
mean free path of the conduction s electrons depends on 
the spin due to s-d scattering and the different d den- 
sity of states at e F for majority and minority spin bands. 
In the present work we consider only the scattering on 
nonmagnetic impurities. 

As a starting point, we take the system of 
Gor'kov equations 18 for the normal and anomalous 
Green's functions G^(xi,x 2 ) = — (T T if;-i (xi)^ (x 2 )) and 
F^(xi,x 2 ) = (T T %l)\(xi)^{x2)), where x = (r, r) is a 



four-component vector and the creation and annihilation 
field operators are associated with s electrons. By carry- 
ing out the Fourier transformation in the xy-plane and 
over the imaginary time r, we get the following system 
for the Green's functions: 
i) for the F-layers: 



1 
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+ e F + h(z) - x -f 2 sd G d J(z, z) 



A*(z)G s Mz,z') + 
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ii) for the S-layer: 
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In Eqs. Q-© K is the in-plane momentum, paral- 
lel to the S/F interface, m is the effective electron mass 
which is assumed to be the same for both metals, h(z) 
is exchange field in the ferromagnet, u> = irT(2n + 1) are 
Matsubara frequencies (the units are % = 1 = &b). The 
scattering processes are introduced in the Born approxi- 
mation. The parameters uq and j sc i are the strengths of 
impurity potentials, c and xo are impurity concentrations 
in the S- and F-layers. We assume that a BCS coupling 
constant is zero for the ferromagnet, therefore A(z) = 
in the F-layers. We also neglect by the possible deviation 
of A(z) from zero in the F-region due to scattering, since 
this correction is of the order of 7 4 d which is small. 

The superconductor order parameter has to be found 
self-consistently, 



A(z) = \tY, 



ndn 
~2^ 



F*(u,k,z = z'), 



(4) 



where summation over uj goes up to Debye frequency ujd , 
A > is the BCS coupling constant in a superconductor, 
and F — F^. The critical temperature T c is defined as 
the first zero of equation A(z) = when T decreases 
from high temperatures. 
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Below in this section and in Sec. Ill we present a 
scheme to evaluate the Green's functions considering as 
the first step the non-self-consistent solution of Eqs. 1)112(1 
where A(z) = A is a real number which does not depend 
on z. Sec. IV is devoted to the self-consistent evaluation 
of A(z). We will assume, that the mutual orientation 
of magnetizations in the F layers is antiparallel (AP), 
therefore h(z) = h > in the Fi-layer, and h(z) = — h 
in the F3-layer. The advantage of the AP configuration 
is that in this case the self-consistency can be achieved 
for real values of A(z) in the S region. The study of the 
influence of the mu tual orientation of magnetizations on 
T c (Refs. Ilil2(12l|) in the framework of the given model 
requires to consider A(z) as a complex valued function. 
These question is beyond the present study and will be 
discussed in the forthcoming publication. However, as 
can be seen further, the general conclusions of the given 
paper are not sensitive to the particular configuration of 
the magnetizations. At the first step we also suppose 
that there is no scattering in the S layer. The scattering 
processes in the S layer (Eq. |3J) are taken into account at 
the last step of the evaluation of the critical temperature 
(Sec. V). 

By introducing the Green's functions 
GH(x 1 ,x 2 ) = -(T-^JOri)^^)) and F^( Xl ,x 2 ) = 
(T T ipi(xi)ipi(x2)) the system of Gor'kov equations can 
be written in the matrix form 26 



ILUl - A 



G F 
F G 



iS(z-z'), 



(5) 



where I is the unit matrix, and A is the (2 x 2)-matrix 
differential operator, the components of which can be 
found by comparing the Eqs. and Eq. JSJ. 

In order to find the matrix Green's function, consider 
the Schrodinger's equation with the Hamiltonian A: 



iwl - A 



*l){z) = 0. 



(6) 



This equation has four linear independent solutions, 

(m=U), 



and 



(p=U)- 



We require that ^v(z) and ip p (z) obey zero boundary 
conditions at the points z = ±a, and choose these in- 
dependent solutions in such a way that two functions 
ipi(z) and ipl( z ) describe spin-up electrons in the ferro- 
magnetic layers, and functions <Pi(z) and ^Pi(z) describe 
spin-down holes in the F-layers. Namely, in the layer Fi 
(—a < z < —d) the solutions f^{z) have the form 



<Pl (z) = I x I sui 



sin 



p[(z + a) 
p{(z + a) 



(7) 



and in the layer F3 [d < z < a) the solutions ip p (z) are 



j sin 


p\{a- 


z) 


I sin 




z) 



(8) 



HerepJ^ is an electron (hole) momentum in the layer Fi, 

(9) 



P\ {1) 



,; 2m ( ep — - — ± h ± — r^Jj ± iu> j . 



2m 



and p 3 are momenta in the layer F3, 



2m \ £F ~2^ Th± 2 T m) ±lw )- (10) 



The inverse life-times of quasiparticles are given by 
r f(j) = -2x o7 2 d ImGff (u) = k F (l) /(ml m) ), here k F il) = 

y/2m{ep ± h) being Fermi momenta in the ferromagnet 
and Zjn) being mean free paths which are considered as 
parameters. 

In the S-region (— d < z < d) the solutions of Eq. © 
are 



e ik_{z+d) _|_ gfl 



-i/c-l- (z-\-d) 



-ik— (z-\-d) 



(ii) 



{z— d) 



+ D% (") e"-'-« + D'_ (") e-> k -<'-'\ 
where the wave vectors k± are defined as 



k± = W2m (e F - — ±iy/u 2 + A 2 
2m 



and 



Q = A 



a/w 2 + A 2 - uj 



We neglect the interfacial roughness, thus the coefficients 
A±, B±, C±, D± have to be found from the conditions 
of continuity of the functions ipp,(z) and ip P (z) and their 
derivatives at the points z — ±d, that can be done easily 
by solving the system of algebraic linear equations. 

To evaluate the matrix Green's function, let us intro- 
duce the matrices 
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and let J be the matrix of " currents" , 



_ ( Jn 3u 

hi hi 



with components 



G(z,z') = 2m $(z)[J- l ] Ti &(z') 



(13) 



Here T denotes the transposition operation. The ob- 
tained expression allows to evaluate the normal and 
anomalous Green's functions in both layers (S and F). 



III. ANOMALOUS GREEN'S FUNCTION 
A. S-layer 

Consider first the anomalous Green's function (Cooper 
pair wave function) F(u>, k, z = z') in the S-region. De- 
note 9± = 2ik±d, and 6± = 9 ± ifi, i.e., 9 and S are 
real and imaginary parts of phases 9± . Using solutions 
of Eq. © in the superconductor, we get the exact 
expressions for currents 



j w = (l~a 2 )j 



where 



f w = 2ik+ [AtC p + e- l9 + -A^Cle 1 ^} 
- 2ik_ [BtD p + e- i6 - 



(14) 



Since the currents j MP do not depend on z, the same 
expressions can be obtained using the solutions of Eq. © 
in the ferromagnetic layers. 

It is convenient to introduce an energy variable £ = 
ep — n 2 /2m. The typical dependence of F(£) on £ un- 
der given arguments cu and A at the point z — z' = 
is shown in Fig. 1. Function F(£) exhibits the quantum 
oscillations which are the result of exponentials e ±i6± in 
Eq. (|14f> with rapidly varying phases. Since the super- 
conducting order parameter is determined by the integral 
of over £, one can average F(£) over the oscillations. 

Denote a pp (/i, p =f , J,) the components of the matrix 



[J 



_ 1]T _ j u TT a n 



For Ou P we get 



a 



sign(/ip) 
pp (l-a 2 )Dcn 



7° 

J — U.. — C 



1 

Dot 



a np c 



AP configuration, m = 
E ex = 0, k p = 0.826 A 1 
"= 1 20 A, I = 40 A 

up down 

A = 9.25 K 



3w = Vz %{z) - <p„ (z) V 2 </> (z), (12) g 



here fi, p =f , J., and Vz= (V z — Vz) is the antisymmetric 
gradient operator. The matrix J is the Wronskian of the 
system (jHJl which does not depend on z, i.e., dJ(z)/dz = 
0. Finally, the matrix Green's function introduced in (JSJ 
is given by 
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FIG. 1: The typical dependence of the anomalous Green's 
function F(lj,£, z) on the energy £ = s 2 /2m- ef under given 
A = 9.25 K, lj = 0, 2 = and e cx = 0, mean free paths in 
the F-layer are l-\ = 120 A and ly = 40 A, Fermi momentum 

h, F = 0.826 A" 1 . F(£) exhibits quantum oscillations, the 
smooth line shows function F(£) averaged over oscillations, 

i. e. Fb(£) (see the text for the details). 



where 



"HP 



sign(^p) 
{l-o?) 

sign(^p) 



2ik+AZ p C~ p e 



-p n~Po—s 



1-a 2 



2ik+Al p CZ p e- s 



(15) 



2ik_B~ p DZ p e s 



here 9 ± id = 9± = 2ik±d, and Den = det J/ (1 — a 2 ) 2 is 
the determinant of the matrix of currents: 



Den 



-D + r+e 



2i8 



r_e' 



-US 



The expressions for Do and T± are given in Appendix A. 

By carrying out the Fourier transformation of a MP , we 
can write the first terms of the expansion: 



b + e 1 ' 



b PP e ~ 



(16) 



where b^ p are defined by the following integrals 



1 

2^ 



r+e 1 * + r_e"^ - D ' 



Using Eq. I|13jl. we get expression for F(ui, £, z = z 1 ): 
F{u,Z,z = z') = 2m 

where ipZ(z) and ip£{z) are the components of solutions 
<Pi_i(z) and ipp(z) in the S-layer. Denote 

#i ± i<$i = k±(d + z), #2 i ^2 = fc±(d — z), 
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then 9i + 9 2 = 0, Si + S 2 = 6 and 
where 

= aA+e~ Sl + B p e Sl , A p _ = aA p e Sl + B p e~ Sl , 
Y, p + = C p _e- &2 + aD p _e 52 , Y, p _ = C p + e &2 + aD p + e' &2 . 

The higher order terms with e , where n > 2, can 
be dropped in the expansion Qltjjl because they are re- 
sponsible for rapid oscillations of F(£). Finally, we come 
to the following expression for the anomalous Green's 
function in the S-layer, 



F{u, z = z') = F (u>, C, z) + Fx(u>, £, z), 



where 



2m y £{KK^ + ^-Etb+ p }, (17) 



pp 



bt„e 2i6 



b PP e 



-2l0 2 



lip 



A P Z P + [, 



pp 



The contribution Fi to the function F is essential only 
in the vicinity of S/F interfaces, z = ±d, as far as 6i ~ 
(z + d) and 02 ~ (d — z). At the point z — z' = (the 
middle of the S-layer) the anomalous Green's function is 
determined by the function Fq(u>, £, z) which is shown by 
the thick smooth line in Fig. 1. The obtained result is 
used below in Sec. IV where we discuss the self-consistent 
evaluation of the order parameter. 



B. F-layer 

Due to proximity effect, the correlations between elec- 
trons are induced in the ferromagnet close to the super- 
coducting layer. Instead of simple decay, as it would 
be for the superconductor/normal-metal interface, in the 
case of ferromagnetic layer the Cooper pair wave func- 
tion exhibits the damping oscillatory behavior in the 
ferromagnet with increasing a distance from the S/F 
interfac e) 3 ! 4 ! 15 . The reason is that exchange splitting of 
bands in the F-region changes the pairing conditions for 
electrons, therefore the Cooper pairs are formed from 
quasiparticles with equal energies but with different in 
modulus momenta and —p\ . Due to the non-zero cen- 
ter of mass momentum Ap, the Cooper pair wave func- 
tion obtains the spatially dependent phase in the ferro- 
magnetic layer. In the "clean" limit (no scattering in the 
ferromagnet) one can fincU^ that the Cooper pair wave 
function oscillates with the distance z into the F-layer as 
~ sm(z/£%)/(z/Q) where = v F /e cx . 

This result holds also in the case of " dirty" ferromag- 
net. The microscopic theory of S/F multilayers based on 



the quasiclassical Usadel equations 3 ^ predicts that the 
anomalous Green's function behaves in the ferromagnet 
as ~ exp{— (1 + i)\fh) Dmz}, where Dm = vpl/3 is the 
diffusion coefficient and I is the electron mean free path 
in the F-layer. Therefore, a length scale for oscillations 
and damping is the same and this scale is set by the 
length £m = \/27£p/3. Below in this section it is shown 
that in the framework of our model the scales for oscil- 
lations and damping of the anomalous Green's function 
arc determined by different lengths. 

We can find the anomalous Green's function 
F(ui, £, z = z') in the F-region (d < z < a) following the 
same approach that was used to evaluate the F-function 
in the superconductor. The solutions tp^n\(z) in the layer 
F3 (d < z < a) are given by Eq. (JHJ). For the solutions 
<fn(z) (p =|, 4) we can write 



-ipUa—z) 



+ Y+ (^^J e ip ^ a_z ) +Y P e -ipi( a - z )^ 

where and Y± can be found from conditions of con- 
tinuity of the functions ip^{z) and their derivatives at 
z = d, assuming perfect S/F interface. 

The anomalous Green's function averaged over oscilla- 
tions is 

F(cj,£,z = z') = 2m^^(z)(a MT )^+(z) 
p 

= 2m^F> MT >. (18) 

p 

It turns out that function f~(z) contains four terms with 
multipliers e ±e + and e ±e - . Denoting 9± — 6±i6, we can 
write F^ in the form 



F, = ^ 



$-e-*° 



(19) 



where 



p 



6+ cos^zi) + sin(^zi) 

x sm\pl(d F - Zx)}, 

= e p cos(j^zi) + E p sin(p^zi) 

x sin[p|(djr - zi)}, 



(20) 



here z\ = z — d is the distance from the S/F interface, 
and 
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e p = aA p _e s ^- n p °- & 



B p _e 

= — B+e 6 +a^±A':, 
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l t 1 
Pi 



Pi 



a !h+A p e s + ^-B p e- s 
v Ps Pi > 
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Cooper pair wave function 



±iw in Eq. llOj l. since |o;| < wp, we obtain 
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FIG. 2: The behavior of the anomalous Green's function F(z) 
along the F/S/F structure with ds = 400 A, d F = 50 A, 
£ cx = 1.156 eV, Z T = 133 A, l x = 35 A. The point z = 200 A 
is the S/F interface. Solid line — Re F(z), dashed dotted 
line — Im F(z). The contribution Re F (z) [see Eq. |TT7J] to 
the function Re F(z) in the S-layer is shown by dashed line. 



Using Eqs. ((T^|) . ljT%|l . ifT^ll we get the expression for 
function F averaged over the rapid oscillations 



F(u,£,z = z') = 2m 



EE 



TO 



It follows from Eq. (|2U|) for that the dependence 
of function F(u>, £, z) on variable z or z\ — z — d is given 
by a sum of the terms with sine and cosine from argu- 
ments (pl+Ps)Zi and [p\—p^)zi. The terms with phases 
(f>3 +p\i) z i determine the short-periodic oscillations with 
respect to oscillations with larger period ~ (pj — pg)" 1 . 
Neglecting the nonessential terms with short-periodic os- 
cillations, the anomalous Green's function can be pre- 
sented in the form 

F(u,£,z = z') = e(w,e)cos(Ap 3 ^i) (21) 
+ S(w,^)sin(Ap 3 zi), 

t i t( n 
where zi = z — d, Ap% = p 3 — pg, and p 3 are given by 

Eqs. ©, P|. 

The real and imaginary parts of the function 

F(z) = T^ I" dZF^z) 

normalized on the value of its real part at the point z = 
(the middle of the S-layer) are shown in Fig. 2. The 
dashed line in Fig. 2 shows the contribution Fo(z) [see 
Eq. p7[l] to the function F(z) in the S-layer. We can 
estimate the lengths responsible for the oscillations and 
decay using Eq. i|l(J|) for momenta p\ andp 3 . Neglecting 



Pl il) = 



yj2m{£ ± h) 



l±-r-i 



4 ^H±h 



As far as the integration over £ goes from till ep, then 
the damping of oscillations is determined by the value 



1 

To 



1 



1 

w 



Neglecting ±i7j,j\ and ±iuj in Eq. (|ll)fl . we get 



= V2K 



ft 



If £ ~ £p, then 



Ap 3 zi 



-zi 



7TZl 
vm 



where ^ = ttvf/e C x — is the half-period of oscil- 
lations of F(z) (the distance between the nearest zeros). 
Note, that ^ irlo in contrast to what is found by using 
a quasiclassical approach. 

The oscillatory terms in Eq. (|21fl arise due to quantum 
interference between two plane waves describing an elec- 
tron and a hole propagating in the ferromagnetic layer 

T I 

with different momenta p 3 and — p% along the z-axis. If 
h ^ then Ap 3 =^ 0, and the oscillatory dependence 
of the Cooper pair wave function occurs due to the ex- 
change field in the ferromagnet. If h = 0, then F(z) 
exhibits only the exponential decay into the F-layer with 
characteristic length l . As it was already pointed out 
by many authors, the physical picture of the proximity 
effect is similar to the nonuniform Fulde-Ferrel-Larkin- 
Ovchinnikov (FFLO) stat o 27 i 28 which is characterized by 
the oscillatory dependent order parameter and arises in a 
homogeneous superconductor in the presence of a strong 
enough uniform exchange field. 



IV. SELF-CONSISTENT EVALUATION OF THE 
ORDER PARAMETER 

In this section we proceed further by constructing the 
self-consistent solution of Eqs. (HIL'I 111 . In case of antipar- 
allel orientation of magnetizations in the ferromagnetic 
layers the self-consistency can be reached if the order pa- 
rameter A(z) takes the real values. We will search the 
self-consistent solution of Eq. (0J in the S-layer assuming 
that in this equation the function F(lj, £, z) is replaced by 
its first contribution Fo(u>,£,z) given by Eq. 1|17(1 . Func- 
tion Fq(z) is shown by dashed line in Fig. 2 and can be 
approximated by a simple analytical function on z like 
oc cos(qz), where q is a parameter. In order to take into 
account the correction F\(z) one has to choose a more 
complex class of sample functions for A(z). However, 
this will not change the results significantly. 
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Let us look for A(z) in the form 
A(z) = Acos(qz) ~ A ( 1 



2 2 

q z 



where the wave vector q (which has to be found) is small. 
The magnitude A(d) = (1— 5q)A defines the amplitude of 
the superconducting order parameter at the S/F interface 
(see Fig. 2), here So — q 2 d 2 /2 ~ 0.1 is a small param- 
eter. Following the well-known WKB-approximation 29 , 
we search the solutions of Schrodinger's equation (01 in 
the form 



i>{z) = 



e '£+(*0 



r2\ 



For £± we get the system of equations 

1 

2m 

1 

2m 



A*{z)e 



*+ + A(z)e*- = 0, (22) 

■■£-£) 



= o, 



where £ = ep — k 2 /2m and primes above £± denote the 
derivatives by z. 

In case of q = 0, A(z) = A, four solutions of sys- 
tem (I22|) are £9 = ±k+z, — ±k+z + iloga and 
£9 = ±fc_z + iloga, £_ = ±/c_z which give four ini- 
tial eigenfunctions of the non-perturbed Eq. © (q = 0): 

1' 



4(z) = 



4M = Q 



Consider, for example, the perturbed solution u + (z) 
which corresponds to vP + {z) in the case of q ^ 0. We 
look for the phases £± in the form £± = £± + ?7±, where 
£9 = k + z, = k + z + iloga for vP + (z). The typical 
order of rj± is ~ qz ~ = V2<5o < 1- By linearizing the 
system (1221) with respect to 7?± we come to the following 
equations 



1 k 4 



V+ + iA{r-/- - r) + ) 



9 2 

q z 



(23) 



a — rf_ — iA(r;_ — 77+) = A 
m 



2 2 

<7 z 



We also have dropped the terms with r]± and r]± which 
are small as compared to /c±t?±, since 77^ ~ V± 1 l±/d ~ 
■q' ± y/28^/d < fc±r4 if d ~ 200 A and fc± - 0.5 A -1 , and 
»7± ~ ?7±/d <C k±rj±. The equations similar to Eqs. f2l^l 
can be written also for phases rj± which determine other 
three solutions u—(z) and v±{z). Solving these equations 
we get 



u±(z) 



with 



-±,3 



Z° + T^Z 2 + T^Z 



T±Z* 



.±„2 



'0 > 



'0 > 



and 



w ± (z) = 



' oe ±i[k-*-K£°(*)r 
o ± l [fc_ z +ci ±) ( z )] 



where 



pfz 3 + P2 z2 + Pi z + Pen 
pf^ 3 _ + Pi z ~ Pa- 



The expressions for coefficients T i and of polynomials 
are given in Appendix B. 

Next the procedure of evaluation of the anomalous 
Green's function in the S-layer is similar to one described 
in details in Sees. II and III for the case of A(z) = A. By 
representing the solutions <fip,(z) and °f ^q. © as 

a linear combination of eigenfunctions u±(z) and v±(z) 
similar to representation Ijlip. we can find the new coef- 
ficients A±, B±, C±, D± solving the system of 4 linear 
equations. By evaluating the currents j^p at the point 
z = (jp P do not depend on z) , we obtain the expressions 
for jup similar to Eq. (|14() where k± should be replaced 
by 



k + = fcj 



1-a 2 



Pi 



a 2 pt 



and 6± = 2ik±d = 9 ± id. The substitutions k± — > k± 
also have to be made in Eq. (|15|l for a^ p and in the expres- 
sion for detJ (see Appendix A). Finally, the anomalous 
Green's function F(u>,£,z) is given by Eqs. l(T7|) where 
A± and are replaced by new functions A± and S^.: 



K 






At 


= aA^ 1 - 1 ^ 


+ B p e- Sl - lC - } 


K 


_ (JP e -&2-iv i f ) 






= C p + e S ' +tr >+ ) 4 


-aD p e- 5 * + < +) 



here S±, 82, JJ± and are functions on z. The fixed 
point q = q* which determines the order parameter A(z) 
has to be found numerically by solving Eq. (0J using the 
iterative procedure. 



V. CRITICAL TEMPERATURE T c 

If the anomalous Green's function in the S-region 
F(oj, £,z) — F(cj,^,z = 0)cos(g*z) is known, the super- 
conducting transition temperature T c can be found. Up 
to now we assume the "clean" limit for a superconductor. 
The corrections (0 due to scattering will be taking into 
account further. Let us introduce the function 

= I * <% F(oj,Z,z = 0), 
kf Jo 
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where kp is Fermi momentum in the S-layer. This in- 
tegral can be evaluated only numerically. However, we 
can approximate by the analytical function of argu- 
ment uj. Let us represent in the form 



F,„ 



VU 1 



For the bulk superconductor F^ = 1. Let T — > T c , 
therefore A <C A(0), where A(0) is the order parameter 
at T = 0. If uj takes values from till 
be well approximated by the following function 



5cl>£>, Fu can 



FP ~ A,tanh 



f 7oM\ 
V 2^d / ' 



(24) 



The coefficients Aq and 70 are found numerically by 
minimizing the norm of a difference between the exact 
and approximate function. These coefficients are non- 
monotonic functions of the F-layer thickness cLf when ds 
is fixed. For typical values of the parameters describing 
the F/S/F structure the magnitudes of A and 70 are 
A ~ 0.9 and 7o ~ 4.0. 

The scattering in the S-layer is introduced by Eq. |J3}. 
Numerical analysis shows that in Eq. © the Green's 
function G^(z, z) does not depend on z in the S-region 
and its real part is negligibly small. Obviously, A w (z) = 
A w cos(<7*z). From numerical analysis it follows that G u 
in the S-layer can be represented as 

= i^ EF ^Gf T ( W ,e,^ = o) 

« -A - 



'OT' (25) 

Taking into account Eqs. I|24l) and l|25l) . equations © 
can be written in the form similar to the case of bulk 
superconductor—: 



OJ « UJ - 



Aq uj' 

2t y^TAj' 



2r v/^ + Aj' 



(26) 



where t^ 1 = 2wcUqN(ef) is the inverse life-time of quasi- 
particles in the superconductor, and N(ep) = mkp/2'K 2 
is density of states at the Fermi energy. Deriving Eq. Ij26|l 
we took into account that, if \t^ 1 ~ ujp, ~ 300 K cor- 
responding to mean free path l s ~ 130 A, then for to' ~ 

+ j4 /2to and 70 ~ 4.0 we have tanh ( r y \uj'\/2uju) » 1. 

Equations (|26|) can be written as2& 



u/ = wr7(w), A w = Arj(u)), 



2t Vuj 2 + A 2 
Using H27|l and (J2J we come to the equation for T c 

J_ tanhf 70 '"' ' 



2ujj 



(27) 



(28) 



where 

/ ^0 
uj=uj + — , p = p ^o<Po, 
2t 

and po = ^N(ep) is the renormalized coupling constant. 
By carrying out the summation over Matsubara frequen- 
cies uj = irT c (2n + 1) in Eq. (J2EJ, we get the equation for 
reduced critical temperature r = T c /T®: 



r = exp 



r. 



%(t) 



(29) 



where 



E 



4r e 



-(2n+l)7j 



£ (2n + l)(l + r e-( 2 «+ 1 )"o) ! 



7 7rT 

%(t) = T, 



T = exp 



2t uj d J ' 



and = 27r 1 a;D7e 



-i/po 



(7 



C = 0.577. 



transition temperature of the bulk superconductor. 



VI. RESULTS AND DISCUSSION 

In this section we present the results of numerical cal- 
culation of the critical temperature T c . We first focus on 
the general features of a behavior of the system. Next 
we consider selected experimental data which can be in- 
terpreted in the framework of the given model. 



A. Oscillatory behavior of T c 

The typical dependence of critical temperature T c (dp) 
with respect to ferromagnetic layer thickness dp with 
ds = 400 A is shown in Fig. 3 where the model parame- 
ters are given in the figure caption. The effective electron 
mass is m = m e (m e is a bare electron mass) . For super- 
conductor we took ojd — 276 K and T® = 9.25 K which 
are the parameters of bulk Nb. The corresponding nor- 
malized magnitude (1 — Sq) of the order parameter at the 
S/F interface as a functions of dp is shown in Fig. 4. 

Both functions (1 — 5q) and T c (dp) show the pro- 
nounced damping oscillatory behavior with the same pe- 
riod. The oscillatory behavior of T c (dp) is a consequence 
of oscillations of the amplitude A(d) = A(l — do) of the 
order parameter at the S/F interface when dp \s vary- 
ing. The minima of A(d) correspond to minima of T c 
and the maxima of A(d) correspond to the maxima of 
T c , as they should. The oscillations of A(d) in turn are 
caused by the oscillations of the anomalous Green's func- 
tion F(z) in the F-layer. Function F(z) must satisfy the 
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FIG. 3: Critical temperature T c for the F/S/F trilayer as a 
function of the ferromagnetic layer thickness dp. The param- 
eters are e ex = 0.771 eV (exchange field in the ferromagnet), 
Z T = 129 A, l± = 37 A, l 3 = 132 A(mean free paths in the 
F and S layers), Uf = 0.826 A" 1 (Fermi momentum in the 
superconductor), ujd = 275 K, T c ° = 9.25 K. 



zero boundary condition at the ferromagnet/ vacuum in- 
terface. Because of oscillations of F(z) in the F-region, 
the order parameter at the S /F interface is forced to ad- 
just in such a way that the condition F(a) = is fulfilled 
at the outer boundary z = a of the F-layer. 

The results of numerical analysis, presented in Table I 
for different values of exchange field e ex and effective elec- 
tron mass m, show that the period £f of T c -oscillations 
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FIG. 4: The dependence of the normalized amplitude (1 — So) 
of the superconducting order parameter at the S/F interface 
as a function of the F-layer thickness (If for the same param- 
eters as in Fig. 3. 



-l 

F ! 



(30) 



here fc^ is the Fermi momentum in a superconductor. 
The period £f, therefore, does not depend on the elec- 
tron mean free paths in the S- and F-layers. The first 
minimum of T c (cIf) occurs at the thickness £f/2, while 
the location of first maximum is 

As can be seen from Fig. 5, the strong scattering in the 
ferromagnetic layers significantly damps the oscillations 
of T e , but their period remains unchanged for any val- 
ues of the mean free paths l-\ and As it follows from 
the analysis presented in Sec. Ill, the reason of such a 
behavior is that the strong scattering in the F-region af- 
fects only the length Iq of decay of the Cooper pair wave 
function F(z) but not the period ~ £^ of its oscillations. 
The less pronounced are the oscillations of F(z) with re- 
spect to z\ = z — d in case of strong electron scattering, 
the less is the amplitude of oscillations of A(d) and T c 
with respect to the ferromagnetic layer thickness dp- In 
case of extremely strong scattering, the coherent coupling 
which was established due to these oscillations between 
two boundaries of ferromagnetic layer is destroyed and 
thus the oscillations of T c are suppressed completely. 

We also observed that strong scattering in the S layer 
(small mean free path l s ) suppresses the amplitude of T c 
oscillations (look at Fig. 6). The critical temperature is 
higher for smaller values of l s . The reason for it is that 
in the thin superconducting films T c is reduced with re- 
spect to T° due to dimensional effect, and the magnitude 
of T c depends on ds only via the dimensionless thick- 
ness ds/^s, where £s oc VCo's is a coherence length for 
the dirty superconductor, £ is a BCS coherence length. 
Small mean free path / s , therefore, corresponds to large 
value of the effective film thickness d$/£,s- 



TABLE I: The comparison of the periods of oscillations £f 
predicted by formula 1301 with the period obtained from 
numerical analysis for different values of exchange field e ex , 
effective electron mass m and a superconductor layer thick- 
ness ds- Other model parameters are the same as in Fig. 3. 
The accuracy of determining of £f is restricted by a finite step 
for dp in numerical calculations. 



Eex (eV) 


m (m e ) 


ds (A) 


(A) 


£f (A) 


0.385 


1.0 


400 


13.97 


14.0 


0.771 


1.0 


400 


9.98 


10.0 


1.156 


1.0 


400 


8.06 
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2.027 


1.0 
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FIG. 5: The critical temperature T c (<1f) for the case of 
weak and strong scattering in the ferromagnetic layers; e cx = 
0.385 eV, L = 132 A. Dots (weak scattering): If = 124 A, 
k = 38 A; squares (strong scattering): If = 32.5 A, Zj. = 9 A. 



B. Comparison with experiment 

Experimental situation on the oscillatory behavior of 
T c (cZf) in the S/F structures is known to be controver- 
sial. Nevertheless, there are two groups of experiments 
described in the literature where oscillations of T c (dp) 
were clearly observed and the 3d fcrromagnets were used 
as F layers — these are reports on Fe/Nb/Fe trilayers by 
Miihge et al&, and Nb/Co and V/Co multilayers by Obi 
et al~. 

In Fig. 7 the fitting is shown to experimental data 
by Muhge et al. 8 for Fe(d F )/Nb(400 A)/Fe(cZ F ) trilay- 
ers prepared by rf sputtering. According to formula 130(1 



E ex = 0.385 eV 

I 6 = 1 24 A, I, = 38 A 

up down 

—m— l s = 132 A 

— o— | S = 265 A 



10 



20 



30 



40 



50 



FIG. 6: The critical temperature T" c (dj?) for different values of 
the mean free paths in the S-layer; e cx = 0.385 eV, i-j = 124 A, 
Zj = 38 A. Squares: l s = 132 A; dots: l s = 265 A. 



the period £ F of oscillations is determined by exchange 
splitting energy in the ferromagnet. If we take the value 
e d cx ~ 0.149 Ry = 2.03 eV (Ref. |p) of exchange split- 
ting of the Fe d-bands near the Fermi energy and put 
to = to 6 , we obtain £^ = 6.09 A (see Table I) which is 
too small as compared to the location of a maximum at 
dp ~ 10 -i- 15 A in Fig. 7. However, we can assume that 
in the S and F layers the Cooper pairs are formed by s 
electrons of Nb and Fe. The value of exchange splitting 
e s cx = 0.028 Ry = 0.381 eV at the bottom of Fe s bands 
(r point, Ref. |3(] ) together with to = m e gives the pe- 
riod £ F = 14.05 A. Thus, the first minimum of T c (dp) 
is at the point £ F /2 rs 7 A, and the first maximum is 
at £ F rs 14 A. From Fig. 7 it follows that these values 
correlate with positions of minimum and maximum of 
T c which can be roughly determined from the scattered 
experimental points. We have put ep = 0.387 Ry corre- 
sponding to the s band of Nb (Ref. I30T) which gives the 
Fermi momentum value kp = 1.18 A -1 for m = m e . We 
used uj d = 276 K and T c ° = 9.25 K for Nb. The fitting 
parameters are the values of mean free paths in Fe and 
Nb which were estimated approximately as l-\ — 120 A, 
l\ = 40 A, and l s = 269 A. Note, that magnetic measure- 
ments by Miihge et al. showed that thin Fe layers were 
not magnetic for dp < 7 A, and it was assumed that mag- 
netically " dead" Fe-Nb alloy of a thickness about 7 A was 
formed at the interfacial S/F region for all samples with 
different d F - Miihge et al. qualitatively explained the 
observed non-monotonic behavior of T c (dp) in terms of 
a rather complex behavior of this magnetically "dead" 
Fe-Nb layer when dp was varying (see details in Ref. @). 
They also argued that a non-monotonic T c (dp) behavior 
in their case could not be possible due to the mechanism 
of 7r coupling as it was predicted for the S/F multilayers 
because of a single S layer in the trilayer system. In- 
deed, the well-known theoretical prediction by Buzdin et 
o/A4 ascribes the oscillatory behavior of T c (dp) to the 
periodical switching of the ground state energy between 
0- and 7r-phases of the order parameter if the neighbor- 
ing S-layers in the S/F multilayer are coupled. However, 
it follows from the above analysis that the oscillatory 
behavior of T c (dp) does not necessary require the 7r cou- 
pling and can occur also for a trilayer (or bilayer) F/S/F 
structure. 

Let us consider the experiments on Nb/Co multilay- 
ers by Obi et al. 9 . The theoretical curve T c (dp) in 
comparsion with experimental data is shown in Fig. 8. 
The exchange splitting of Co spin-up and down s bands 
at r point is e cx = 0.014 Ry (Ref. HJ which gives 
= 19.87 A ( to = to g ). The first and second min- 
imum of T c {dp) should, therefore, be placed at points 
^min = 10 A and G^ in = 30 A. These values correlate 
with values 12 A and 32 A obtained from experiment. 
The fitting mean free paths are 1^ = 240 A, l± = 80 A, 
and l s = 188 A. We have to note that in experiment 
Nb/Co structures are multilayers. A qualitative resem- 
blance of theoretical T c curve calculated for a trilayer 
structure with experimental points for a multilayer and 
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FIG. 7: The comparison of the theoretical T c (cIf) curve with 
experiment by Miihge et al. (Ref. H ) for Fe/Nb(400A)/Fe 
trilayers. The fitting parameters are if = 120 A, l± = 40 A, 
l s = 269 A. 



the agreement between theoretical and experimental val- 
ues of d^ in and d^ in allows us to assume that neighboring 
S layers were decoupled in the experiment. As it was ob- 
served by Strunk et alJSL for similar Nb/Fe multilayered 
system (where F-layer is 3c? transition metal), the decou- 
pling regime is set when dp is larger than some critical 
value g?c° which in turn is less than the critical thick- 
ness df of the onset of ferromagnetism. This threshold 
value was w 7 A in experiments by Obi et al& In 
Ref. Oj it was noted that df was less than the first min- 
imum of T c at w 12 A, so that for Nb/Co system 
the first minimum could not be ascribed to the onset of 
ferromagnetism as it was argued by Miihge et al. for the 
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FIG. 8: The comparison of the theoretical T c (cIf) curve with 
experiment by Obi et al. (Ref. 0) for Nb(400A)/Co multi- 
layers. The fitting parameters are Z T = 240 A, ij, = 80 A, 
l s = 188 A. 



Fe/Nb/Fe system^. Our theoretical explanation assum- 
ing the decoupling regime is incorrect only for very thin 
Fe layers with dp < dj? c when, probably, the Fe films are 
nonmagnetic due to alloying effect. 

Note also, that experiments by Obi at el£ on 
Nbi-^Ti^/Co multilayers with Nbi-^Ti^ alloy being su- 
perconductor with small coherence length did not reveal 
the oscillatory behavior of T c but showed only a small 
reduction of the critical temperature T c w 8 K for large 
dp as compared to the bulk value T° w 9.2 K (see Fig. 3 
in Ref. 0). Therefore, the observation of increasing of 
T c when the scattering is strong in S-layer together with 
damping of oscillations for small l s (see Fig. 6) is in a 
qualitative agreement with these experimental observa- 
tions. 



VII. SUMMARY 

In conclusion, we have presented a theory of proxim- 
ity effect in F/S/F trilayer nanostructures where S is a 
superconductor and F are layers of 3d transition ferro- 
magnetic metal. As a starting point of our calculations, 
we took the system of Gor'kov equations, which deter- 
mine the normal and anomalous Green's functions. The 
solution of these equations was found together with a 
self-consistent evaluation of the superconductor order pa- 
rameter. In accordance with the known quasiclassical 
theories of proximity effect for S/F multilayers2^ii^i&, 
we found that due to a presence of exchange field in 
the ferromagnet the anomalous Green's function F(z) ex- 
hibits damping oscillations in the F-layer as a function 
of a distance z from the S/F interface. In the presented 
model a half-period of oscillations of F(z) is determined 
by the length = irvp/e C x, where vf is the Fermi ve- 
locity, e cx is the exchange field, and the length of damp- 
ing is given by la = (1/if + 1/Zj.) -1 , where if and l± 
are spin-dependent mean free paths in the ferromagnetic 
layer. The oscillations of the anomalous Green's func- 
tion (Cooper pair wave function) in the F-region and 
a zero boundary condition at the ferromagnet /vacuum 
interface give rise to the oscillatory dependence of the 
superconductor order parameter at the S/F interface vs 
the F-layer thickness dp. These oscillations result in os- 
cillations of the superconductor transition temperature 
T c (dp) with a period £p = ir/y/me cx . Thus we have 
demonstrated that the nonmonotonic oscillatory depen- 
dence of critical temperature T c (dp) does not necessarily 
require the mechanism of 7r-coupling between neighbor- 
ing superconducting layers as it takes place in the S/F 
multilayers^. The strong electron scattering either in 
the superconductor or in the ferromagnet significantly 
suppresses the oscillations. In case of extremely strong 
scattering in the ferromagnet the length of damping Iq 
becomes very short and the oscillations of T c are sup- 
pressed completely. The reason of that is the loss of 
coherent "coupling" between two boundaries of ferro- 
magnetic layer that was established due to oscillations of 
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Cooper pair wave function F[z). We compared our re- 
sults with existing data on T c (dp) for Fe/Nb/Fe trilayers^ 
and V/Co multilayers^, where F's are 3d ferromagnets, 
and found reasonable agreement with theory and exper- 
iment. 



In case of q ^ four linear independent solutions of 
Eq. |[BJ| have the form: 
i) solution u+iz): 



u + (z) = 
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APPENDIX A 

The determinant Den = det J/(l — a 2 ) 2 of the matrix 
of currents (Eq. 112(1 is given by the expression 



Den = -D Q + T+e 
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and A±, _B±, C±, D± are coefficients introduced in 
Eq. (TTJ. 



where 



r] { +\z) = - iX+W+z 3 + iX+V+z 2 
+ 2i\ 2 + V+z + 2iX\V+ 

= T+Z 3 +T+Z 2 +T+Z + T, 



„«(*) 



3 a, -r , 2 * -r 1 1 * t / , 

1 (+),■. A ( m 
V+ 0) + 



3~ Z 3 + T 2 z2 +T^Z + Tq ; 



q 2 z 3 



ii) solution U-(z): 
U-{z) 

where 



ae 



-ik+z—irk (z) 



iii) solution v+(z) 



= T 3 z3 ~t±z 2 + t 1 ± z-t ± ; 



where 



/ ae ik-z+i( ( + } (z) \ 
V +(z) = l k-Z+i < l _+\z) I 



C| +) (z) = iA_iy_z 3 + ?;a_f_z 2 
- 2iA 2 _y_z + 2U 3 _V r _ 
= p+z 3 + p+z 2 + p+z + p+, 

= p~z 3 + p^z 2 + p7z + p~; 



APPENDIX B 

Let us define the quantities 



A ± = — 



W± = 



v± = 



2 J. a 2 
A 2 / m \ 2 



m 



Vuj 2 + A 2 



iv) solution i>_(z): 
u_(z) 

where 



s — ik-z-iC' + ' (z) ' 
,-ifc_z-iC£~ 3 (2) 



cr j (z) = -ci +) (-^) 



4- o + 2 ± ± 

p 3 z - p 2 z + p x z - p . 
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